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Abstract 

We show that the quantum query complexity of detecting if an n-vertex graph contains a triangle 
is 0{'n?^'^). This improves the previous best algorithm of Belovs [2] making 0(n^^^^^) queries. For the 
problem of determining if an operation o : 5* x 5 — >■ 5* is associative, we give an algorithm making 
QQ^rjio/T-j qygj-jgg^ tjjg fjj-g^ improvement to the trivial 0(|5'|'^''^) application of Grover search. 

Our algorithms are designed using the learning graph framework of Belovs. We give a family of 
algorithms for detecting constant-sized subgraphs, which can possibly be directed and colored. These 
algorithms are designed in a simple high-level language; our main theorem shows how this high-level 
language can be compiled as a learning graph and gives the resulting complexity. 

The key idea to our improvements is to allow more freedom in the parameters of the database kept 
by the algorithm. As in our previous work 9:, the edge slots maintained in the database are specified 
by a graph whose edges are the union of regular bipartite graphs, the overall structure of which mimics 
that of the graph of the certificate. By allowing these bipartite graphs to be unbalanced and of variable 
degree we obtain better algorithms. 

1 Introduction 

Quantum query complexity is a black-box model of quantum computation, where the resource measured is 
the number of queries to the input needed to compute a function. This model captures the great algorithmic 
successes of quantum computing like the search algorithm of Grover [S] and the period finding subroutine of 
Shor's factoring algorithm |12j . while at the same time is simple enough that one can often show tight lower 
bounds. 

Recently, there have been very exciting developments in quantum query complexity. Reichardt 
showed that the general adversary bound, formerly just a lower bound technique for quantum query com- 
plexity 0, is also an upper bound. This characterization opens a new avenue for designing quantum query 
algorithms. The general adversary bound can be written as a relatively simple semidefinite program, thus 
by providing a feasible solution to the minimization form of this program one can upper bound quantum 
query complexity. 

This plan turns out to be quite difficult to implement as the minimization form of the adversary bound 
has exponentially many constraints. Even for simple functions it can be challenging to directly write down 
a feasible solution, much less worry about finding a solution with good objective value. 

To surmount this problem, Belovs [2] introduced the beautiful model of learning graphs, which can be 
viewed as the minimization form of the general adversary bound with additional structure imposed on the 
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form of the solution. This additional structure makes learning graphs easier to reason about by ensuring 
that the constraints are automatically satisfied, leaving one to worry about optimizing the objective value. 

Learning graphs have already proven their worth, with Belovs using this model to give an algorithm 
for triangle finding with complexity 0{nP^/^^), improving the quantum walk algorithm |10| of complexity 
0{n^'^). Belovs' algorithm was generalized to detecting constant-sized subgraphs [131 1^, giving an algorithm 
of complexity o(n^~^^'') for determining if a graph contains a fc- vertex subgraph H, again improving the |10| 
bound of 0(n^~^/'^). All these algorithms use the most basic model of learning graphs, that we also use in 
this paper. A more general model of learning graphs (introduced, though not used in Belovs' original paper) 
was used to give an o(n'^/^) algorithm for fc-element distinctness, when the inputs are promised to be of a 
certain form [3] . Recently, Belovs further generalized the learning graph model and removed this promise to 
obtain an o{rr^^'^) algorithm for the general fc-distinctness problem [1]. 

In this paper, we continue to show the power of the learning graph model. We give an algorithm for 
detecting a triangle in a graph making C>(n^/^) queries. This lowers the exponent of Belovs algorithm from 
about 1.296 to under 1.286. For the problem of testing if an operation o : S x S ^ S is associative, where 
\S\ = n, we give an algorithm making 0{n^*^/'^) queries, the first improvement over the trivial application 
of Grover search making 0(n^/^) queries. Previously, Dorn and Thierauf |3] gave a quantum walk based 
algorithm to test if o : 5 x 5 — >■ 5" is associative that improved on Grover search but only when |S"| < n^^^. 

More generally, we give a family of algorithms for detecting constant-sized subgraphs, which can possibly 
be directed and colored. Algorithms in this family can be designed using a simple high-level language. Our 
main theorem shows how to compile this language as a learning graph, and gives the resulting complexity. 
We now explain in more detail how our algorithms improve over previous work. 

Our contribution. We will explain the new ideas in our algorithm using triangle detection as an 
example. We first review the quantum walk algorithm of |10j , and the learning graph algorithm of Belovs [21 . 
For this high-level overview we just focus on the database of edge slots of the input graph G that is maintained 
by the algorithm. A quantum walk algorithm explicitly maintains such a database, and the nodes of a learning 
graph are labeled by sets of queries which we will similarly interpret as the database of the algorithm. 

In the quantum walk algorithm [TU] the database consists of an r-element subset of the n-vertices of G 
and all the edge slots among these r-vertices. That is, the presence or absence of an edge in G among a 
complete r-element subgraph is maintained by the database. In the learning graph algorithm of Belovs, the 
database consists of a random subgraph with edge density < s < 1 of a complete r-element subgraph. In 
this way, on average, O(sr^) many edge slots are queried among the r-element subset, making it cheaper to 
set up this database. This saving is what results in the improvement of Belovs' algorithm. Both algorithms 
finish by using search plus graph collision to locate a vertex that is connected to the endpoints of an edge 
present in the database, forming a triangle. 

Zhu [l3j and Lee et al. |9] extended the triangle finding algorithm of Belovs to finding constant sized 
subgraphs. While the algorithm of Zhu again maintains a database of a random subgraph of an r-vertex 
complete graph with edge density s, the algorithm of Lee et al. instead used a more structured database. 
Let H he a A:- vertex subgraph with vertices labeled from [k]. To determine if G contains a copy of H, the 
database of the algorithm consists oi k — 1 sets Ai, . . . , A^-i of size r and for every {i, j'} E H — {k} the 
edge slots of G according to a sr-regular bipartite graph between Ai and Aj. Again both algorithms finish 
by using search plus graph collision to find a vertex connected to edges in the database to form a copy of H. 

In this work, our database is again the edge slots of G queried according according to the union of regular 
bipartite graphs whose overall structure mimics the structure of H. Now, however, we allow optimization 
over all parameters of the database — we allow the size of the set Ai to be a parameter r^ that can be 
independently chosen; similarly, we allow the degree of the bipartite graph between Ai and Aj to be a 
variable dij. This greater freedom in the parameters of the database allows the improvement in triangle 
finding from 0(n'^^/^^) to 0{n^^'^). Instead of an r-vertex graph with edge density s, our algorithm uses as a 
database a complete unbalanced bipartite graph with left hand side of size ri and right hand side of size r2 . 
Taking ri < r2 allows a more efficient distribution of resources over the course of the algorithm. As before, 
the algorithm finishes by using search plus graph collision to find a vertex connected to endpoints of an edge 
in the database. 



2 



The extension to functions of the form / : [q']"^" — > {0, 1}, hke associativity, comes from the fact that the 
basic learning graph model that we use depends only on the structure of a 1-certificate and not on the values 
in a 1-certificate. This property means that an algorithm for detecting a subgraph H can be immediately 
applied to detecting H with specified edge colors in a colored graph. 

If an operation o : Sy.S S is non-associative, then there are elements a, 6, c such that ao(boc) 7^ {aob)oc. 
A certificate consists of the 4 (colored and directed) edges boc = e,aoe,aob = d, and doc such that 
a o e ^ d o c. The graph of this certificate is a 4-path with directed edges, and using our algorithm for this 
graph gives complexity 0(151^°/^). 

We provide a high-level language for designing algorithms within our framework. The algorithm begins 
by choosing size parameters for each Ai and degree parameters for the bipartite graph between Ai and Aj. 
Then one can choose the order in which to load vertices ai and edges {ai,aj) of a 1-certificate, according to 
the rules that both endpoints of an edge must be loaded before the edge, and at the end all edges of the 



certificate must be loaded. Our main theorem Theorem 8 shows how to implement this high-level algorithm 



as a learning graph and gives the resulting complexity. 

With larger subgraphs, optimizing over the set size and degree parameters to obtain an algorithm of 
minimal complexity becomes unwieldy to do by hand. Fortunately, this can be phrased as a linear program 
and we provide code to compute a set of optimal parametertj^ 



2 Preliminaries 

The quantum query complexity of a function /, denoted Q{f), is the number of input queries needed to 
evaluate / with error at most 1/3. We refer the reader to the survey [B] for precise definitions and background. 

For any integer q > 1, let [q] = {1,2,..., q}. We will deal with boolean functions of the form / : [g]"^" ^■ 
{0, 1}, where the input to the function can be thought of as the complete directed graph (possibly with self- 
loops) on vertex set [n], whose edges are colored by elements from [q]. When q — 2, the input is of course just 
a directed graph (again possibly with self- loops). A partial assignment is an element of the set {[q] U j*})"^". 
For partial assignments ai and a2 we say that ai is a restriction of 02 (or alternately 02 is an extension of 
ai) if whenever ai(i, j) 7^ * then ai{i,j) = a2{i,j)- A 1-certificate for / is a partial assignment a such that 
/(x) = 1 for every extension x € [g]"^" of a. If a is a 1-certificate and x G [g]"^" is an extension of a, we 
also say that a is a 1-certificate for / and x. A 1-certificate a is minimal if no proper restriction of a is a 
1-certificate. The index set of a 1-certificate a for / is the set = {(i,j) € [n] x [n] : a{i,j) ^ *}. Besides 
these standard notions, we will also need the notion of the graph of a 1-certificate. For a graph G, let V{G) 
denote the set of vertices, and E{G) the set of edges of G. 

Definition 1 (Certificate graph). Let a be a 1-certificate for f : [q]"^" {0, 1}. The certificate graph Ha 
of a is defined by E{Ha) — la, and V{Ha) is the set of elements in [n] which are adjacent to an edge in la- 
The size of a certificate graph is the cardinality of its edges. A minimal certificate graph for x, such that 
f{x) = 1, is the certificate graph of a minimal 1-certificate for f and x. The 1-certificate complexity of f is 
the size of the biggest minimal certificate graph for some x such that f{x) = 1. 

Intuitively, if x £ [q]"^" is an extension of a 1-certificate a, the certificate graph of a represents queries 
that are sufficient to verify /(x) = 1. 

Vertices of our learning graphs will be labeled by sets of edges coming from the union of a bunch of 
bipartite graphs. We will specify these bipartite graphs by their degree sequences, the number of vertices on 
the left hand side and right hand side of a given degree. The following notation will be useful to do this. 

Definition 2 (Type of bipartite graph). A bipartite graph between two sets Yi and Y2 is of type 
{{{ni,di),...,{nj,dj)},{{mi,gi),...,{nii,ge)}) if Yi has Ui vertices of degree di for i ^ 1, . . . , j , and Y2 
has nii vertices of degree gi ^or i — 1,...,£, and this is a complete listing of vertices in the graph, i.e. 
\Yi\ = Yli=i "^i IF2I = T.i=i "^i- ^o^e fl^'So that Yll=i i^i^i = Z]i=i "m-igi- 

^code is available at https://github.com/troyjlee/learning_graph_lp 
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Learning graphs Wc now formally define a learning graph and its complexity. We first define a learning 
graph in the abstract. 

Definition 3 (Learning graph). A learning graph Q is a 5-tuple {y,£,w,^,{py : y G Y}) where (V,£) is 
a rooted, lueighted and directed acyclic graph, the weight function w : £ maps learning graph edges to 

positive real numbers, the length function £ : £ — > N assigns each edge a natural number, and py : £ ^ R is 
a unit flow whose source is the root, for every y ^ Y . 

A learning graph for a function has additional requirements as follows. 

Definition 4 (Learning graph for a function) . Let f : [g]"^" — > {0,1} be a function. A learning graph Q for 
f is a 5-tuple {V,£,S,w,{py : y G /~^(1)}), where S : V 2"^" maps v e V to a label S{v) C [n] x [n] of 
variable indices, and {V,£,w,i,{py : y G /^^(l)}) "is a learning graph for the length function £ defined as 
i{{u,v) — \S{v) \ S{u)\ for each edge {u,v). For the root r G V we have S{r) = 0, and every learning graph 
edge e = {u,v) satisfies S{u) C S{v). For each input y G f^^{l), the set S{v) contains the index set of a 
1-certificate for y on f, for every sink v £ V of py. 

In our construction of learning graphs we usually define S by more colloquially stating the label of each 
vertex. Note that it can be the case for an edge {u, v) that S{u) — S{v) and the length of the edge is zero. 
In Belovs pi what we define here is called a reduced learning graph, and a learning graph is restricted to 
have all edges of length at most one. 

In this paper we will discuss functions whose inputs are themselves graphs. To prevent confusion we will 
refer to vertices and edges of the learning graph as L-vertices and L-edges respectively. 

We now define the complexity of a learning graph. For the analysis it will be helpful to define the 
complexity not just for the entire learning graph but also for stages of the learning graph Q. By level d of Q 
we refer to the set of vertices at distance d from the root. A stage is the set of edges of G between level i and 
level j, for some i < j. For a subset 1^ C V of the L-vertices let = {{v, w) G £ : v G V} and similarly let 
= {{u,v) G £ : V £ V}. For a vertex v we will write instead of {w}^, and similarly for instead of 

M-. 

Definition 5 (Learning graph complexity). Let Q be a learning graph, and let E Q £ be the edges of a stage. 
The negative complexity of E is 

Co(i?) = ^^(e)u;(e). 

eeE 

The positive complexity of E under the flow py is 
The positive complexity of E is 

Ci{E) = maxCiy(£;). 

yeY 

The complexity of E is C{E) = y/CQ{E)Ci{E), and the learning graph complexity of Q is C{Q) = C{£). 
The learning graph complexity of a function f , denoted CQ{f), is the minimum learning graph complexity of 
a learning graph for f . 

Theorem 1 (Belovs). Q{f) = 0{Cg{f)). 

Originally Belovs showed this theorem with an additional log q factor for functions over an input alphabet 
of size q; this logarithmic factor was removed in [3]. 



4 



Analysis of learning graphs Given a learning graph G, the easiest way to obtain another learning graph 
is to modify the weight function of G- We will often use this reweighting scheme to obtain learning graphs 
with better complexity or complexity that is more convenient to analyze. When G is understood from the 
context, and when w' is the new weight function, for the edges i? C £ of a stage, we denote the complexity 
of E with respect to w' by (E). 

The following useful lemma of Belovs gives an example of the reweighting method. It shows how to upper 
bound the complexity of a learning graph by partitioning it into a constant number of stages and summing 
the complexities of the stages. 

Lemma 2 (Belovs). If £ can be partitioned into a constant number k of stages Ei, . . . , Ek, then there exists 
a weight function w' such that 

C^'{G)^0{CiE,) + ... + CiEk)). 

Now we will focus on evaluating the complexity of a stage. Our learning graph algorithm for triangle 
detection is of a very simple form, where all L-edges present in the graph have weight one, all L-vertices in 
a level have the same degree, incoming and outgoing flows are uniform over a subset of i-vertices in each 
level, and all L-edges between levels are of the same length. In this case the complexity of a stage between 
consecutive levels can be estimated quite simply. 

Lemma 3. Consider a stage of a learning graph between consecutive levels. Let V be the set of L-vertices 
at the beginning of the stage. Suppose that each L-vertex v G V is of degree-d with all outgoing L-edges e of 
weight w{e) — 1 and of length £(e) < £. Furthermore, say that the incoming flow is uniform over L-vertices 
W C V , and is uniformly directed from each L-vertex v € W to g of the d possible neighbors. Then the 

complexity of this stage is at most 

Proof. The total weight is d\V\. The flow through each of the g\W\ many L-edges is (glVFj)^^. Plugging 
these into |Definition 5] gives the lemma. □ 

To analyze the cost of our algorithm for triangle detection, we will repeatedly use [Lemma "3} The 
contributions to the complexity of a stage are naturally broken into three parts: the length the vertex 
ratio \V\/\W\, and the degree ratio d/g. This terminology will be helpful in discussing the complexity of 
stages. 



For our more general framework given in Section 4 flows will no longer be uniform. To evaluate the 



complexity in this case, we will use several lemmas developed in |9] . The main idea is to use the symmetry 
of the function to decompose flows as a convex combination of uniform flows over disjoint edge sets. A 
natural extension of |Lemma"3| can then be used to evaluate the complexity. To state the lemma we first 
need a definition. For a set of L-edges E, we let Py{E) denote the value of the flow py over E^ that is 

Definition 6 (Consistent flows). Let E be a stage of G between two consecutive levels, and let Vi, . . . ,Vs be 
a partition of the L-vertices at the beginning of the stage. We say that {py} is consistent with , . . . , 
if Py{V^'^) is independent of y for each i. 

The next lemma is the main tool for evaluating the complexity of learning graphs in our main theorem, 
[Theorem 8l 

Lemma 4 ([9]). Let E be a stage of G between two consecutive levels. Let V be the set of L-vertices at the 
beginning of the stage and suppose that each v d V has outdegree d and all L-edges e of the .stage satisfy 
w{e) — 1 and £{e) < £. Let Vi, . . . ,Vs be a partition of V, and for all y and i, let Wy^i C Vi be the set of 
vertices in Vi which receive positive flow under py . Suppose that 

1. the flows {py} are consistent with {Vi^}, 

2. \Wy^i\ is independent from y for every i, and for all v € Wy^i we have py{v^) = Py{Vi^)/\Wy^i\, 



5 



3. there is a g such that for each vertex v € Wyj the flow is directed uniformly to g of the d many 
neighbors. 

Then there is a new weight function w' such that 




C"' (E) < max^J-p^^ . (1) 



We will refer to max^ as the maximum vertex ratio. For the most part wc will deal with the 

problem of detecting a (possibly directed and colored) subgraph in an n- vertex graph. We will be interested 
in symmetries induced by permuting the elements of [n], as such permutations do not change the property 
of containing a fixed subgraph. We now state two additional lemmas from [9 that use this symmetry to help 
establish the hypotheses of [Lemma 4| 

For a € Sn, we define and also denote by a the permutation over [n] x [n] such that (7{i,j) ~ (c(j), cr(i)). 
Recall that each L- vertex u is labeled by a fc-partite graph on [n], say with color classes Ai, . . . , A^, and 
that we identify an L-vertex with its label. For u e S'„ we define the action of cr on u as <7{u) = v, where v 
is a /c-partite graph with color classes cr(^i), . . . , a{Ak) and edges {cr(i), cr(j)} for every edge {i, j} in u. 

Define an equivalence class [u] of L- vertices by [u] — {cr(w) : a G Sn}- We say that 5„ acts transitively 
on flows {py} if for every y,y' there is a t G Sn such that py{{u,v)) = ((r(u), t(w)) for all L-edges {u,v). 

The following lemma from ^ shows that if Sn acts transitively on a set of flows {py} then they are 
consistent with [v\^ , where v is a vertex at the beginning of a stage between consecutive levels. This will set 



us up to satisfy hypothesis (1) of Lemma 4 



Lemma 5 (j9 ). Consider a learning graph Q and a set of flows {py} such that Sn acts transitively on {py}. 
Let V he the set of L -vertices of Q at some given level. Then {py} is consistent with {[u]"*" : u S V}, and, 
similarly, {py} is consistent with {[u]~ : u G V}. 



The next lemma gives a sufficient condition for hypothesis (2) of Lemma 4 to be satisfied. The partition 



of vertices in Lemma 4 will be taken according to the equivalence classes [u] 



Lemma 6 (0). Consider a learning graph and a set of flows {pj,} such that Sn acts transitively on {py}. 
Suppose that for every L-vertex u and flow py such that py{u^) > 0, 

1. the flow from u is uniformly directed to g^{[u]) many neighbors, 

2. for every L-vertex w, the number of incoming edges with from [w] to u is g~{[w], [u]). 

Then for every L-vertex u the flow entering [u] is uniformly distributed over Wj, [„] C [u] where [„]| is 
independent ofy. 



3 Triangle algorithm 

Theorem 7. There is a bounded-error quantum query algorithm for detecting if an n-vertex graph contains 
a triangle making 0{n^/'^) many queries. 

Proof. We will show the theorem by giving a learning graph of the claimed complexity, which is sufficient by 
[Theorem 1[ We will define the learning graph by stages; let Vt denote the L-vertices of the learning graph 
present at the beginning of stage t. The L-edges between Vt and Vt+i are defined in the obvious way — there 
is an L-edge between vt £ Vt and Vt+i G Vt+i if the graph labeling vt is a subgraph of the graph labeling 
Ut+i, and all such L-edges have weight one. The root of the learning graph is labeled by the empty graph. 

For a positive input graph G, let ai, 02, 03 be the vertices of a triangle of G. The algorithm (see Figure [T]) 
depends on set size parameters ri,r2 G [n\,ri,r2 — o{n) and a vertex degree parameter A G [n] that will be 
optimized later. We will choose ri < r2 such that r2/ri is an integer. The cost of each stage will be upper 



bounded using Lemma 3 
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Stage 1 (Setup): The initial level Vi consists of the root of the learning graph labeled by the empty 
graph. The level V2 consists of all L-vertices labeled by a complete unbalanced bipartite graph with disjoint 
color classes Ai,A2 C [n] where \Ai\ — ri ~ 1 and \A2\ = r2 — 1 and ri <r2- Flow is uniform from the root 
to all -L-vertices such that ^ Ai, <^ A2 for i = 1, 2. 3. 



Cost: The hypotheses of Lemma 3 hold trivially at this stage. The length of this stage is 0{rir2)- The 



vertex ratio is 1, and the degree ratio is 



("r)r 



-ri- 
r2 



= 0(1), 



as ri,r2 = o(n). Thus the overall cost is 0{rir2)- 



Stage 2 (Load oi): During this stage we add a vertex to the set Ai and connect it to all vertices in ^2- 
Formally, V3 consists of all vertices labeled by a complete bipartite graph between color classes Ai , A2 of 
sizes ri, r2 — 1, respectively. The flow goes uniformly to those L-vertices where ai is the vertex added to Ai. 
Cost: By the definition of stage 1, the flow is uniform over L-vertices at the beginning of stage 2. The 
out-degree of every L-vertex in Vi is n — ri — r2 + 2. Of these, in L-vertices with flow, exactly one edge is 



taken by the flow. Thus we can apply Lemma 3 Since the degree ratio was 0(1) for the flrst stage, the 
vertex ratio is also 0(1) for this stage. The length is r2 — 1. The degree ratio is 0{n). Thus the cost of this 
stage is 0{^/nr2)- 



Stage 3 (Load 02): We add a vertex to A2 and connect it to all of the ri many vertices in Ai. Thus the 
L-vertices at the end of stage 3 consist of all complete bipartite graphs between sets Ai,A2 of sizes ri,r2, 
respectively. The flow goes uniformly to those L-vertices where 02 is added at this stage to ^2- Note that 
since we work with a complete bipartite graph, if ai e Ai and 02 G A2 then the edge {ai, 02} is automatically 
present. 

Cost: The amount of flow in a vertex with flow at the beginning of stage 3 is the same as at the beginning of 
stage 2, as the flow out-degree in stage 2 was one and there was no merging of flow. Thus flow is still uniform 
at the beginning of stage 3. The out-degree of each L-vertex is n — ri — r2 -I- 1 and again for L-vertices with 
flow, the flow out-degree is exactly one. Thus we can again apply [Lemma 3 

The length of this stage is ri. The vertex ratio is 0{n/ri) as flow is present in L-vertices where ai is in 
the set Ai of size ri (and such that 02,03 are not loaded which only affects things by a 0(1) factor). The 
degree ratio is again 0{n) as the flow only uses L-edges where 02 is added out of 71 — ri — r2 + 1 possible 
choices. Thus the cost of this stage is 0{\J nlr\yfnr\) = 0{72y/ri). 

Stage 4 (Load 03): We pick a vertex v and A many edges connecting v to A2. Thus the L-vertices at the 
end of stage 4 are labeled by edges that are the union of two bipartite graphs: a complete bipartite graph 
between Ai,A2 of sizes ri, r2, and a bipartite graph between and A2 of type {(1, A)}, {(A, 1), (r2 — A, 0)}. 
Flow goes uniformly to those L-vertices where u = 03 and the edge {02, 03} is not loaded. 
Cost: Again the amount of flow in a vertex with flow at the beginning of stage 4 is the same as at the 
beginning of stage 3, as the flow out-degree in stage 3 was one and there was no merging of flow. Thus the 
flow is still uniform. The out-degree of L-vertices is {n — ri~ r2)(J'^) , and the flow out-degree is {^^^^) ■ Thus 
we can again apply [Lemma 3[ 

The length of this stage is A. At the beginning of stage 4 flow is present in those L-vertices where 
tti S Ai,a2 G A2 and 03 is not loaded. Thus the vertex ratio is 0{{n/ri){n/r2))- Finally, the degree ratio is 
0(n). Thus the overall cost of this stage is 
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Stage 5 (Load {02,03}): We add one new edge between v and A2. Thus the L-vertices at the end of 
this stage will be labeled by the union of edges in two bipartite graphs: a complete bipartite graph between 
Ai,A2 of sizes ri, r2, and the second between v and A2 of type {(1, A + 1)}, {(A + 1, 1), (r2 — A — 1, 0)}. Flow 
goes uniformly along those L-edges where the edge added is {02, 03}. 

Cost: The flow is uniform at the beginning of this stage, as it was uniform at the beginning of stage 4, the 
flow out-degree was constant in stage 4, and there was no merging of flow. Each L-vertex has out-degree 
r2 — A and the flow-outdegree is one. Thus we can again apply [Lemma 3[ 

The length of this stage is one. The vertex ratio is 0{{n/ri){n/r2)n) as flow is present in a constant 
fraction of those L-vertices where ai G Ai,a2 & A2 and v = 03. The degree ratio is r2 — A, as there are this 
many possible edges to add and the flow uses one. Thus the overall cost of this stage is 




Stage 6 (Load {01,03}): We add one new edge between v and Ai. Thus the L-vertices at the end of 
this stage will be labeled by the union of three bipartite graphs between Ai,A2 and v,A2 as before, and 
additionally between v, Ai of type {(1, 1)}, {(1, 1), (ri — 1, 0)}. Flow goes uniformly on those L-edges where 
{oi, 03} is added. 

Cost: Again flow is uniform as it was at the beginning of stage 5, the flow out-degree was constant and 
there was no merging. Each L-vertex has out degree ri and the flow out-degree is one. Thus we can again 
apply [Lemma "3} 

The length of this stage is one. The vertex ratio is 0{{n/ri){n/r2)n{r2/X)) as flow is present in a constant 
fraction of those L-vertices where oi S Ai,02 G A2,v = and {02,03} is present. The degree ratio is ri. 
Thus the overall cost of this stage is 

By choosing ri = n^/"^ , r2 = n^/^, A = n'^/^ we can make all costs, and thus their sum, 0(n^/^). □ 

To quickly compute the stage costs, it is useful to associate to each stage a local cost and global cost. The 
local cost is the product of the square root of the degree ratio and the length of a stage. The global cost is 
the square root of the factor by which the stage increases the vertex ratio — we call this a global cost as it is 
propagated from one stage to the next. Thus the square root of the vertex ratio at stage t will be given by 
the product of the global costs of stages 1, . . . , i — 1. As the cost of each stage is the product of the square 
root of the vertex ratio, square root of the degree ratio, and length, it can be computed by multiplying the 
local cost of the stage with the product of the global costs of all previous stages. 



Stage 


1 


2 


3 


4 


5 


6 


Global cost 


1 


\J'nlrx 


\/n/r2 


\/n 


v/^2/A 




Local cost 


rir2 


\/nr2 








V^i 


Cost 


rir2 


\/nr2 


riy/ri 


n'''^A/y/rir2 




n3/VVA 


Value 


^9/7 


^17/14 


„9/7 


„9/7 


„17/14 


^9/7 



4 An abstract language for learning graphs 

In this section we develop a high-level language for designing algorithms to detect constant-sized subgraphs, 
and more generally to compute functions / : [q]"^" — {0,1} with constant-sized 1-certificate complexity. 
This high-level language consists of commands like "load a vertex" or "load an edge" that makes the algorithm 
easy to understand. Our main theorem, [Theorem 8[ compiles this high-level language into a learning graph 
and bounds the complexity of the resulting quantum query algorithm. After the theorem is proven, we can 
design quantum query algorithms using only the high-level language, without reference to learning graphs. 
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This saves the algorithm designer from having to make many repetitive arguments as in Section 3 and also 
allows computer search to find the best algorithm within our framework. 



4.1 Special case: subgraph containment 

We now give an overview of our algorithmic framework and its implementation in learning graphs. We first 
use the framework for computing the function fn ■ [2] {0, 1}, which is by definition 1 if the undirected 

n- vertex input graph contains a copy of some fixed fc-vertex graph H = ([k], E{H)) as a subgraph. This 
case contains all the essential ideas; after showing this, it will be easy to generalize the theorem in few more 
steps to any function / : [q]^^) {0, 1} or / : [q]"^" — {0, 1} with constant-sized 1-certificate complexity. 

Fix a positive instance x, and vertices ai,...,afc S [n] constituting a copy of H in x, that is, such 
that x^ai.aj} = 1 8-11 {hj} S E{H). Vertices of the learning graph will be labeled by fc-partite graphs 
with color classes Ai, . . . , Ak- The sets Ai, . . . ,Ak are allowed to overlap. Each L- vertex label will contain 
an undirected bipartite graph G^ijj — (^niin{ij}i ^max{ij}i for every edge {i,j} G E{H), where 

E{i,j} C A^in{z,j} X ^max{ij}- For {ij} S E{H), by {a.i,aj} we mean (a^^aj) if i < j, and {aj^a^) if j < i. 
For an edge G E{H), and u G [n], the degree of u in Gij towards Aj is the number of vertices in Aj 

connected to u ii u d Ai, and is otherwise. The edges of these bipartite graphs define naturally the input 
edges formally required in the definition of the learning graph: for u ^ v, both (u, w) and {v,u) define the 
input edge {u, v}. We will disregard multiple input edges as well as self loops corresponding to edges (u, u). 
Observe that various L- vertex labels may correspond to the same set of input edges. For the ease of notation 
we will denote G^i jy by both Gij and Gji. We will use similar convention for E^i jy which will be denoted 
by both Eij and Eji. 

Our high-level language consists of three types of commands. The first is a setup command. This is 
implemented by choosing sets Ai, . . . , Ak Q [n] of sizes ri, . . . , rj, and bipartite graphs Gij between Ai and 
Aj for all E E{H). Both the set sizes ri,. .. ,rk and the average degree of vertices in the bipartite 

graph between Ai and Aj are parameters of the algorithm. The degree parameter dij — dji represents the 
average degree of vertices in the smaller of Ai , Aj towards the bigger one in Gij . It is defined in this fashion 
so that it is always an integer and at least one — the average degree of the larger of Ai , Aj can be less than 
one. Without loss of generality there is only one setup step and it happens at the beginning of the algorithm. 

The other commands allowed are to load a vertex ai and to load an edge {ai,aj} corresponding to 
G E{H) (this terminology was introduced by Belovs). There are two regimes for loading an edge. One 
is the dense case, where all vertices in the graph Gij have a neighbor; the other is the sparse case, where 
some vertices in the larger of Ai, Aj have no neighbors in the smaller. We need to separate these two cases 
as they apparently have different costs (and cost analyses). The algorithm is defined by a choice of set sizes 
and degree parameters, and a loading schedule giving the order in which the vertices and edges are loaded 
and which loads all edges of H. 

We now define the parameters specifying an algorithm more formally. 

Definition 7 (Admissible parameters). Let H — ([k], E{H)) he a k-vertex graph, ri, . . . , e [n] he set size 
parameters, and dij € [n] for G E{H) he degree parameters. Then {ri}, {dij} are admissible for H if 

• ^ <ri < n/4 for all i e [k], 

• 1 < dij < mBx{ri,rj} for all {i,j} G E{H), 

• for all i there exists j such that {i,j} G E{H) and dij{2rj + \)/{2ri + 1) > 1. 

We give a brief explanation of the purpose of each of these conditions. We will encounter terms of the 
form (".)/("~'^) that we wish to be 0(1); this is ensured by the first condition. As dij represents the average 
degree of the vertices in the smaller of Ai,Aj towards the larger, the second condition states that this degree 
cannot be larger than the number of distinct possible neighbors. The third item ensures that the average 
degree of vertices in Ai is at least one in the bipartite graph with some Aj . 
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Definition 8 (Loading scliedule). Let H = {[k], E{H)) be a k-vertex graph with m edges. A loading schedule 
for H is a sequence S = S1S2 ■ ■ ■ st+m whose elements Si G [k] or Si E E(H) are vertex labels or edge labels 
of H such that an edge {«, j} only appears in S after i and j, and S contains all edges of H. Let VSt be the 
set of vertices in S before position t and similarly ESt the set of edges in S before position t. 

We can now state the main tlieorem of tliis section. 

Theorem 8. Let H = {[k], E{H)) be a k-vertex graph. Let ri, . . . ,rk,dij be admissible parameters for H, 
and S be a loading schedule for H . Then the quantum query complexity of determining if an n-vertex graph 
contains H as a subgraph is at most a constant times the maximum of the following quantities: 

• Setup cost: 

min{ru,ry}duv, 

{u,v}eE{H) 



Cost of loading St ~ i: 




ma.x{ru,r^} 



X 



.]}eE(H) j:{ij}eE(H) 

ri<rj ri>rj 

Cost of loading St = {i,j} in the dense case where (2niin{ri,rj} + > (2 max{ri, Tj} + 1); 



rjd,j 




max{ri, rj}, 



Cost of loading St = {i,j} in the sparse case where (2min{ri,rj} + l)c?ij < (2 maxjri, rj} + 1); 




max{?\^} 

dill! 



If {i,j} is loaded in the dense case we call it a type 1 edge, and if it loaded in the sparse case we call it 
a type 2 edge. The costs of a stage given by |Thcorem 8| can again be understood more simply in terms of 
local costs and global costs. We give the local and global cost for each stage in the table below. 



Stage 


Global Cost 


Local Cost 




Setup 


1 


E{„.^,}ei^Inin{'■«,^ 




Load vertex i 




y/n X total degree 


of i 


Load a type 1 edge {i,j} 


^max{ri,rj}/dij 


max{ri,rj} 




Load a type 2 edge {i,j} 


^max{ri,rj}/dij 







Proof. We show the theorem by giving a learning graph of the stated complexity. Vertices of the learning 
graph will be labeled by fc-partite graphs with color classes Ai, . . . , oi cardinality (of order) ri, . . . ,rk G [n]. 
The parameter dij > 1 is the average degree of vertices in the smaller of Ai,Aj towards the bigger in the 
bipartite graph Gij. 

The bipartite graph Gij, for each edge {i,j} G E{H), will be specified by its type, that is by its degree 
sequences as given in [Definition 2[ 

We first need to modify the set size parameters {ri} to satisfy a technical condition. Let r;^ < • ■ ■ < ri^, 
be a listing in increasing order. We set r^^ — ri^ and r[^ = 0(rij such that (2r^^ + l)/(2r^^_^ + 1) is an odd 
integer. As a consequence, (2 max{r^, r^} + l)/(2 min{r,-, r'j\ + 1) is an odd integer, for every i 7^ j. We now 
suppose this is done and drop the primes. 

Throughout the construction of the learning graph we will deal with two cases for the bipartite graph 
between Ai and Aj, depending on the size and degree parameters. 
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Figure 2: Example of a part of learning graph corresponding to Case 1 and restricted to the bipartite graph 
between Ai and Aj, where ri < rj. Observe that A,j w -^dij. The loading schedule is 'setup', 'load i', 'load 
j' and 'load 
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A-i : 2ri vertices 




degree 1 



1 new vertex 




Ai with 1 new vertex — >■ 2ri + 1 vertices 



Aj : 2rj vertices 




Ai : 2ri + 1 vertices 
all of degree dij 




Aj with 1 new vertex — )- 2rj + 1 vertices 



Ti vertex-disjoint new edges 




Ai : 2ri + 1 vertices 
all of degree dij 




2rjl+ 1 - (2ri + l)dij v*tices 
with dec ree 




Aj : 2rj + 1 vertices 



Ai : 2ri + 1 vertices 



Aj : 2rj + 1 vertices 



Figure 3: Similar to Figure [2j but for Case 2. 
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Case 1 is where (2 min{ri , } + l)dij > 2max{ri,rj} + 1, which means that there are enough edges 
from the smaller of Ai, Aj to cover the larger. We will say that the parameters for are of type 1. 

In this case, we take d-^ = ®{dij) to be such that 



for some integer Xij. This can be done as (2 maxfr^, r-,} + l)/(2min{ri, rj} + 1) is an odd integer. 
In our construction, Xij will be the average degree of the vertices in the larger of Ai,Aj towards the 
smaller, which we want to be integer. We now consider this done and drop the primes. 

• Case 2 is where (2min{ri,rj} + l)dij < 2max{ri,rj} + 1. We will say that the parameters for 
are of type 2. In this case, all degrees of vertices in the larger of Ai,Aj towards the smaller will be 
either zero or one. 

Now we are ready to describe the learning graph. Figures [2] and |3] illustrate the evolution of a learning 
graph for a subsequence of some loading schedule, that is the sequence of instructions 'setup', 

'load i\ 'load j' and 'load The figures only represent the added edges between Ai and Aj, where 

ri < Tj. Figure [2] corresponds to Case 1, and Figure [3] to Case 2. 

Recall that for every positive instance x, we fixed ai,...,ak G [n] be such that x^a^^a^j ~ 1 for all 
{u,v} € E{H). During the construction we will specify for every edge {u,v} G E{H), and for every stage 
number t, the correct degree cd{u,v,t) which is the degree of in towards a„ in each L-vertex of Vt+i 
with positive flow. 



Stage (Setup): For each edge £ E{H) we setup a bipartite graph between Ai and Aj. The type of 

the bipartite graph depends on the type of the parameters for {i, j}. Let t = min{rj, r^} and g = maxfr^, rj}. 



Case 1: Solving for in Equation (2) we get — {{2£ + l)dij + £ — g)/{2g + 1). Intuitively, dij 
represents the average degree of vertices in the smaller of Ai , Aj and Xij the average degree in the 
larger. Formally, the type of bipartite graph between Ai , Aj , with the listing of degrees for the smaller 
set given first, is ({(2^ - Xij,dij), {Xij, dij - l)},{{2g - d,j,Xij), {dij, Xij - 1)}). 

Case 2: In this case the type of bipartite graph between Ai and Aj, with the listing of degrees for the 
smaller set given first, is {{{2£, d^j)},{{2£dij, 1), {2g - 2£d^j,0)}). 



The L- vertices at the end of stage will be labeled by (possibly overlapping) sets Ai, . . . ,Ak of sizes ri, . . . , 
and edges corresponding to a graph of the appropriate type between Ai and Aj for all G E{H). Flow 

goes uniformly to those L- vertices where none of ai,...,ak are in any of the sets Ai, . . . , A^- For all 
{u, v} e E{H), we set cd(u, v, 0) — 0. 



Stage t when St — i: In this stage we load a^. The L-edges in this stage select a vertex v and add it to 
Ai. For all j such that {i, j} G E{H) we add the following edges: 

• Case 1: Say the parameters for are of type 1. If < rj, then v is connected to those vertices of 
degree Xij — 1 in Aj, and we set cd{i,j,t) — dij. Otherwise v is connected to those vertices of degree 
dij — 1 in Aj, and we set cd{i,j,t) = Xij. 

• Case 2: Say the parameters for are of type 2. If < rj then v is connected to dij vertices of 
degree in Aj, and we set cd{i,j,t) = dij. Else no edges are added between v and Aj, and we set 
cd{i,j,t) = 0. 

For all other {u, v), we set cd(u, v, t) = cd{u, v,t—l). Flow goes uniformly on those L-edges where v — ai. 
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Stage t when st = {hj}'- In this stage we load {ai,aj}. Again we break down according to the type of 
the parameters for Let £ — mm{ri,rj} and g = niax{rj,rj}. 

• Case 1: As both and aj have been loaded, between Ai and Aj there is a bipartite graph of type 
{{{2£ + 1, dij)}, {{2g + 1, Xij)}), with the degree listing of the smaller set coming first. If we simply 
added {ui^aj} at this step, and Uj would be uniquely identifiable by their degree and blow up the 
complexity of later stages. 

To combat this, loading {a^, aj} will consist of two substages t.I and t.II. The first substage is a hiding 
step, done to reduce the complexity of having {a^, aj} loaded. Then we actually load {a^, aj}. 

Substage t.I: Let h — {2g + l)/{2£ + 1). We select £ vertices in the smaller of Ai,Aj, and to each of 
these add h many neighbors. All neighbors chosen in this stage are distinct. Thus at the end of this 
stage the type of bipartite graph between Ai and Aj is {{{£, dij + h), {£ + 1, dij)}, {{£{2g + l)/(2£ + 
1), \ij + 1), {{2g + 1)(1 — £/{2£ + 1)), Xij)}). Flow goes uniformly along those L-edges where neither ai 
nor aj receive any new edges. For all {m, v} e E(H), we set cd(u, v,t + 1) ^ cd{u, v, t). 

Substage t.II: The L-edges in this substage select a vertex u in the smaller of Ai, Aj of degree dij and 
add h many neighbors of degree Xij. Flow goes uniformly along those L-edges where u S {0^,0^} and 
{ai,aj} is one of the edges added. Let s be the index of the smaller of the sets Ai,Aj, and let b the 
other index. We set cd(s, b,t + 1) = dij + h, cd(6, s,t + 1) = Xij + 1 and cd(u, v,t + 1) = cd{u, v, t) for 
{u,v} ^ {i,j}. 

• Case 2: As both and Oj have been loaded, there is a bipartite graph of type {{{2£+ l,dij)},{{{2£ + 
l)dij, 1), {2g + 1 — (2^ + l)dij,0)}). We again first do a hiding step, and then add the edge {ai,aj}. 

Substage t.I: We select £ vertices in the smaller of Ai, Aj and to each add a single edge to a vertex of 
degree zero in the larger of Ai, Aj. Flow goes uniformly along those L-edges where no edges adjacent 
to ai,aj are added. For all {u,v} £ E{II), we set cd{u,v,t+ 1) — cd{u,v,t). 

Substage t.II: A single edge is added between a vertex in the smaller of Ai,Aj of degree dij and 
a vertex in the larger of Ai,Aj of degree zero. Flow goes along those L-edges where {ai,aj} is 
added. Let again s be the index of the smaller of the sets Ai,Aj, and let h the other index. We 
set cd(s, h,t+l) = dij + 1, cd(6, s,t + 1) = 1 and cA{u, v,t+l) = cd(w, v, t) for {it, v} ^ {», j}- 

This completes the description of the learning graph. 



Complexity analysis We will use Lemma 4 to evaluate the complexity of each stage. First we need to 



establish the hypothesis of this lemma, which we will do using [Lemma's] and Lemma "6] Remember that 



given CT € Sn, we defined and denoted by cr the permutation over [n] x [n] such that <j(i,f) = {a{i), a{j)). 
First of all let us observe that every ct S is in the automorphism group of the function we are computing, 
since it maps a 1-certificate into a 1-certificate. As the fiow only depends on the 1-certificate graph, this 
implies that Sn acts transitively on the flows and therefore we obtain the conclusion of |Lcmma 5[ 

Let Vt stand for the L-vertices at the beginning of stage t. For a positive input x, and for an L- 
vertex P € Vt, we will denote the incoming flow to P on a; by Px{P) and the number of outgoing edges 
from P with positive flow on x by g^{P). For an L- vertex R £ Vt-i we will denote by g^ fi{P) number 
of incoming edges to P from L-vertices of the isomorphism type of R with positive flow on x, that is 
g^ r{P) = !{■'' G Sn : Px{t{R),P) ^ 0}|. The crucial features of our learning graph construction are the 
following: at every stage, for every L- vertex P and every u e S'„, the L- vertex (j{P) is also present. The 
outgoing flow from an L-vertex is always uniformly distributed among the edges getting flow. The flow 
depends only on the vertices in the input containing a copy of the graph H, and therefore the values gt{P) 
and g~ fi{P), for Px{P) non-zero, depend only on the isomorphism types of P and R. Mathematically, this 
last property translates to: for all t, for all P G Vt, for all R E Vt-i, for all positive inputs x and y, for all 
(7 e 5'n, we have 

[p,(F)^Oandp,(a(P))^0] =^ [gt {P) = g+ {a{P)) g^^iP) ^ glj,{a{P))]. (3) 
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which is exactly the hypothesis of Lemma 6 



Now we have estabhshed the hypotheses of Lemma 4 and turn to evaluating the bound given there. The 



main task is evaluating the maximum vertex ratio of each stage. The general way we will do this is to 
consider an arbitrary vertex P of a stage. We then lower bound the probability that a{P) is in the flow for 
a positive input x and a random permutation a d Sm without using any particulars of P. This will then 
upper bound the maximum vertex ratio. We use the notation P to denote that L- vertex P has at least 
one incoming edge with flow on input x. 

Lemma 9 (Maximum vertex ratio). For any L-vertex P G Vt+i and any positive input x 

pr[a(P) G F.] = I n ^ n 



max{r„, r^} 

Proof. We claim that an L-vertex P in Vt+i, that is at the end of stage t, has flow if and only if 

Vi G VS(, y{i,j} G ESf, we have Ui G Ai and {ai,aj} G Eij, (4) 
Vi G [fc] \ VSt, V{i, j} G H{E) \ ESt, we have a, ^ Ai and {a^, a,} ^ Eij, (5) 
j} G ESt, the degree of in towards Aj is cd{i,j,t). (6) 

The only if part of the claim is obvious by the construction of the learning graph. The if part can be proven 
by induction on t. For t = 0, the first half ^ is exactly the one which defines the flow for L- vertices in Vi. 

For the inductive step let us suppose first that St ~ i. Consider the label P' by dropping the vertex 
from Ai. Then in P' every bipartite graph is of appropriate type for level t because of ([6]), and therefore 
P' G Vt. It is easy to check that P' also satisfies all three conditions, (for (|6| we also have to use the second 
half of ([5]): {ai,aj} ^ Eij), and therefore has positive fiow. Since P' is a predecessor of P is the learning 
graph, P has also positive flow. 

Now let us suppose that St — {i,i}. In P the edge set Eij can be decomposed into the disjoint union 
of El U E2, where Ei a bipartite graph of type ({(2£ + l,dij)},{{2g + l,Aij)}) and E2 is of type ({(£ + 
(^,0},{((^ + l)/i, 1), {2g + l~{l+l)h,Q)}), and ^ implies that {a,,aj} G E2. Consider the label P' by 
dropping the edges of E2 from Eij. Again, P' satisfies the inductive hypotheses, and therefore gets positive 
fiow, which implies the same for P. 

Suppose now that the P- vertex P is labeled by sets Ai, . . . ,Ak (some may be empty) and let the set of 
edges between Ai and Aj be Eij. We want to lower bound the probability that cr(P) G F^., meaning that 
(t(P) satisfies the above three conditions. Item ([5]) is always satisfied with constant probability; moreover, 
conditioned on item ([s]) the probability of the other events does not decrease. Thus we take this constant 
factor loss and focus on the items Q, ([6]). 

We also claim that, conditioned on item Q holding, item ([6| holds with constant probability. This can 
be seen as in the hiding step, in both case 1 and case 2, the probability that a,;,aj have the correct degree 
given that they are loaded is at least 1/4. In the step of loading an edge, again in case 1 half the vertices 
on the left and right hand sides have the correct degree and so this probability is again 1/4; in case 2, given 
that the edge is loaded, whichever of ai,aj is in the larger set will automatically have the correct degree, 
and the other one will have correct degree with probability 1/2. Now we take this constant factor loss to 
obtain that Prg- [cr(P) G F^] is lower bounded by a constant factor times the probability that item Q holds. 

The events in the first condition are independent, except that for the edge {0^,0^} to be loaded the 
vertices ai and Oj have to be also loaded. Thus we can lower bound the probability it is satisfied by 

Pr[a(P) G P,] - rif TT Pr[a, G (t(A,)] X TT Pr[{a„ a,} G a(P„)|a, G a(A;), a, G C7(A,)" 

ievSt (u,v)eESt 

Now Pro-[ai G f7(Ai)] — Q{ri/n) as this fraction of permutations will put ai into a set of size r^. For the 
edges we use the following lemma. 
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Lemma 10. LetYi,Y2C[n] be of size £, g respectively, and let {yi,y2) E Yi XY2. Let K be a bipartite graph 
betweenYi andY2 of type {{£,d)},{{g,£d/g)}. Then PTc[{yi,y2} E cr{K) — d/ g. 

Proof. Because of symmetry, this probability does not depend on the choice of {2/1,2/2}; denote it by p. Let 
Ki, Kc he a.n enumeration of all bipartite graphs isomorphic to K. We will count in two different ways 
the cardinality x of the set {(e, /i) : e e A';,}. Every contains £d edges, therefore x = cid- On the other 
hand, every edge appears in pc graphs, therefore x = (-gpc, and thus p = d/g. □ 

In our case, the graph Gij as in the hypothesis of the lemma plus some additional edges. By monotonicity, 
it follows that 

Pr[{aj,aj} e a{Eij)\ai e a{Ai),aj e u{Aj)\ = fl{dij/ max{rj,rj}). 

(7 

□ 

This analysis is common to all the stages. Now we go through each type of stage in turn to evaluate the 
stage specific length and degree ratio. 



Setup Cost: The length of this stage is upper bounded by 

inm{ri,rj}dij. 

We can upper bound the degree ratio by 

nS<2'==o(i) 

ie[k] V 2ri ) 

as Tj < n/4. 



Stage t when St = i: In a stage loading a vertex the degree ratio is 0{n) as there are n — ri possible 
vertices to add yet only one is used by the flow. The length of this stage is the total degree which is upper 
bounded by 

ri<rj ri>rj 

Stage t when .sj = {i,j}'. Technically we should analyze the complexity of the two substages as two 
distinct stages. However, as we will see, in both cases the degree ratio in the first substage is 0(1), and 
therefore the local cost of this stage is just the maximum of the local cost of the two substages. 



Stage t.I: In Case 1, the length of this stage is 0(max{ri, rj}) and the degree ratio is constant. In Case 2, 
the length of this stage is 0(min{r,,rj}) and the degree ratio is constant. 



Stage t.IT. In Case 1, the length is /i = 0(max{ri, ^j}/ minjri, rj}). The degree ratio is of order £ 




0{£'^). Thus the square root of the degree ratio times the length is of order max{ri,rj}. 

In Case 2, the length is one and the degree ratio is 0{rirj) as there are 0{rirj) many possible edges that 
could be added and the flow uses one. 

Thus in Case 1 in both substages the product of the length and square root of degree ratio is 
0(max{rj, Tj}). In Case 2, substage II dominates the complexity where the product of the length and 
square root of degree ratio is 0{y/fifj). □ 
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4.2 Extensions and basic properties 



We now extend Theorem 8 to the general case of computing a function / : [q]"'^" {0, 1} with constant- 
sized 1-certificates. A certificate graph for such a function will be a directed graph possibly with self-loops. 
Between i and j there can be bidirectional edges, that is both and (j, i) present in the certificate graph, 
but there will not be multiple edges between i and j, as there are no repetitions of indices in a certificate. 



We start off by modifying the algorithm of Theorem 8 to work for detecting directed graphs with possible 



self-loops. To do this, the following transformation will be useful. 

Definition 9. Let H be a directed graph, possibly with self-loops. The undirected version U{H) of H is a 
simple undirected graph formed by eliminating any self-loops in H , and making all edges of H undirected and 
single. 

Lemma 11. Let H he a directed k-vertex graph, possibly with self loops. Then the quantum query complexity 
of detecting if an n-vertex directed graph G contains H as a subgraph is at most a constant times the 



complexity given in Theorem 8 of detecting U{H) in an n-vertex undirected graph. 



Proof. Let H he a. directed fc- vertex graph (possibly with self- loops) and H' = U{H) be its undirected 
version. Let ri, . . . , rk,dij be admissible parameters for H' , and S a loading schedule for H' . Fix a directed 
n-vertex graph G containing H as a subgraph. Let oi, . . . , be vertices of G such that (a^, aj) G E{G) for 



(i, j) € E{H). We convert the algorithm for loading H' in Theorem 8 into one for loading H of the same 
complexity. 

The setup step for H' is modified as follows. In the bipartite graph between Ai and Aj, if both 
[h G E{H) then all edges between Ai and Aj are directed in both directions; otherwise, if 

{i,j) € E{H) or {j,i) E E(H) they are directed from Ai to Aj or vice versa, respectively. For every 
self-loop in H, say {i,i) G E{H), we add self-loops to the vertices in Ai. Note that these modifications at 
most double the number of edges added, and hence the cost, of the setup step. 

Loading a vertex: When loading Oi we connect it as before, now orienting the edges according to 
or (j, i) in E{H), or both. If {i, i) S E{H), then we add a self loop to a^. The only change in the complexity 
of this stage is again the length, which at most doubles. Notice that in the case of a self-loop we have also 
already loaded the edge {ai,ai). We do not incur an extra cost for loading this edge, however, as the self 
loop is loaded if and only if the vertex is. 

Loading an edge: Say that we are at the stage where St — € E{H'). If exactly one of [i,]), {j, i) G 

E{H) then this step happens exactly as before, except that the bipartite graph has edges directed from 
Ai to Aj or vice versa, respectively. If both and {j,i) G E{H), then in this step all edges added are 
bidirectional. This again at most doubles the length, and does not affect the degree flow probability as 
(oi, Oj) is loaded if and only if {oj, ai) is loaded as all edges are bidirectional. □ 

Lemma 12. Let f : [g]"^" — > {0, 1} be a function such that all minimal 1-certificate graphs are isomorphic 
to a directed k-vertex graph H . Then the quantum query complexity of computing f is at most the complexity 
of detecting H in an n-vertex graph, as given by \Lemma 11\ 

Proof. We will show the theorem by giving a learning graph algorithm. Let Q = {y,£,S,w,{py\) be the 



learning graph from Lemma 11 for H. All oiV,£,S,w will remain the same in our learning graph Q' for /. 
We now describe the definition of the flows in Q' . 

Consider a positive input x to f, and let a be a minimal 1-certificate for x such that the certificate graph 
Ha is isomorphic to H. The flow p^ will be defined as the flow for (thought of as an n-vertex graph, thus 
with n — k isolated vertices) in Q, the learning graph for detecting H. This latter flow has the property that 
the label of every terminal of flow contains E{Ha) and thus will also contain the index set of a 1-certificate 
for X. 

The positive complexity of the learning graph for / will be the same as that for detecting H and the 
negative complexity will be at most that as in the learning graph for detecting H, thus we conclude that the 
complexity of computing / is at most that for detecting H as given in [Lemma 11[ □ 
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Theorem 13. Say that the 1-certificate complexity of f : [g]"^" — > {0, 1} is at most a constant m, and 
let Hi , . . . , He be the set of graphs ( on at most m edges ) for which there is some positive input x such that 
Hi is a minimal 1-certificate graph for x. Then the quantum query complexity of computing f is at most a 
constant times the maximum of the complexities of detecting Hi for i — 1, . . . , c as given by \Lemma 11\ 

Proof. Consider learning graphs Qi, . . . ,Qc given by |Lemina 11 [ for detecting . . . , H^ respectively. Further 
suppose these learning graphs are normalized such that their negative and positive complexities are equal. 

We construct a learning graph Q for / where the edges and vertices are given by connecting a new root 
node by an edge of weight one to the root nodes of each oi Qi, . . . ,Qc- Thus the negative complexity of Q is 
at most c(l + maxi C^^Qi). 

Now we construct the flow for a positive input x. Let a be a minimal 1-certificatc for x such that the 
certificate graph is isomorphic to Hi, for some i. Then the flow on x is first directed entirely to the 
root node of Qi. It is then defined within Qi as in [Lemma 12[ Thus the positive complexity of Q is at most 

c(l + max,Ci(5,))- n 



To make [Theorem 8 and Lemma 11 easier to apply, here we establish some basic intuitive properties 



about the complexity of the algorithm for different subgraphs. Namely, we show that if H' is a subgraph 



of H then the complexity given by Lemma 11 for detecting H' is at most that of H. We show a similar 
statement when H' is a vertex contraction of H. 

Lemma 14. Let H be a directed k-vertex graph (possibly with self-loops) and H' a subgraph of H . Then the 
quantum query complexity of determining if an n-vertex graph G contains H' is at most that of determining 



if G contains H from Lemma 11 



Proof. Assume that the vertices of H are labeled from [k] and that H' is labeled such G E{H) for all 
{i,3)€E{H'). 

The learning graph we use for detecting H' is the same as that for H. For a graph G containing a H' as 
a subgraph, let ai, . . . ,afc be such that {ai,aj) G G for all (i, j) G H' . (If t is an isolated vertex in if', then 
at can be chosen arbitrarily). The flow for G is defined in the same way as in the learning graph for H . Note 
that once ai, . . . , afe have been identified, the definition of flow depends only edge slots — not on edges — thus 
this deflnition remains valid for H' . Furthermore all terminals of flow are labeled by edge slots {ai,aj) for 
all {i,j) G H, and so also contain the edge slots for H' . Thus this is a valid flow for detecting H' . As the 
learning graph and flow are the same, the complexity will be as that given in [Lemma 11[ □ 

Lemma 15. Let H be a k-vertex graph and H' a vertex contraction of H . Then the quantum query complexity 



of detecting H' is at most that of detecting H given in Lemma T7| 

Proof. Again we assume that the vertices of H are labeled from [k] . The key point is the following: if H' is a 
vertex contraction of H, then there are zi, . . . , Zk G [k] (not necessarily distinct) such that {zi, zj) G E{H') if 
and only if {i,j) G E{H). The learning graph for H' will be the same as that for H except for the flows. For 
a graph G containing H' , we choose vertices ai, . . . , (not necessarily distinct) such that if (z^, Zj) G E{H') 



then (ai, Oj) G E{G). As (zi, Zj) G E(H') if and only if {i,j) G E{H), we can define the flow as in Lemma 11 
for ai, . . . , Cfe to load a copy of H' . (Note that there is no restriction in the proof of that theorem that the 
sets Ai,. . . ,Ak be distinct). This gives an algorithm for detecting H' with complexity at most that given 
by [Lemma ll[ for detecting H. □ 



5 Associativity testing 

Consider an operation o : S" x 5 — > S* and let n = \S\. We wish to determine if o is associative on S, meaning 
that ao (bo c) = (aob) o c for all a,b,c € S. We are given black box access to o, that is, we can make queries 
of the form (a, b) and receive the answer ao b. 

Theorem 16. Let S be a set of size n and o : S x S ^ S be an operation that can be accessed in black-box 
fashion by queries (a, b) returning a o b. There is a bounded-error quantum query algorithm to determine if 
{o,S) is associative making 0{n^^/'^) queries. 
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''°£ ,ao(boc) ° gob hoc (a o b) o c ° ° 

Certificate <J=> a o (b o c) {a o b) o c 



= oi =05 
Certificate -i=i> (02 o 03 = ai, 03 o 04 = 05 and 02 o ai 7^ as o 04) 



Figure 4: The 5-vertex certificate graph for associativity. Both pictures represent the same graph certificate 
H, where the second one has been labehed according to the notations of our abstract language. 



Proof. If o is not associative, then there is a triple 02,03,04 such that 02 o (03 o 04) 7^ (02 o 03) o 04. A 
certificate to the non- associativity of o is given by 03 o 04 = oi, 02 o oi, 02 o 03 = 05, and 05 o 04 such that 
02 o oi ^ 05 o 04 (see Figure |4|. Note that not all of oi, . . . , 05 need to be distinct. 

Let H hea directed graph on 5 vertices with directed edges (2, 1), (2, 3), (3, 4), (5, 4). Each non-associative 
input has a certificate graph either isomorphic to or a vertex contraction of H, in the case that not all of 
oi, . . . , 05 are distinct. By [Lemma 15 the complexity of a detecting a vertex contraction of H is dominated 
by that of detecting H, and so by Theorem 13|it suffi ces to show the theorem for H. 



We use the algorithmic framework of 



,6/7 



,5/7 



, r5 = 1 and ^2 



26/7 d 



Theorem 
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to load the graph H. Let ri 
dsi = ri^/^,(i54 = 1. Here di 



n/10,r2 — n*/^,r3 = 
indicates the average 



degree of vertices in the smaller of Ai,Aj for edges directed from Ai to Aj. It can be checked that this 
is an admissible set of parameters. Note that as r^d^^ << r^, loading 05 o 04 will be done in the sparse 
regime. We use the loading schedule S = [1, 2, 4, 3, (2, 1), (2, 3), (3, 4), 5, (5, 4)]. The setup cost becomes 
'^2^21 + ''2^23 + ''4^34 + ^^5^54 = n^"/'', and the costs of loading the vertices and edges are all bounded by 
n^^/'^ as given in the following tables. 



load Oi 
\/nlri 
/nr2d2i/ri 



load 02 
_\/n/r~2 
y/n{d2i + ^23) 



load 04 

^/n(d34 + rsd^i/rj) 
n3/2 , , 7— 



load 03 

Vn{r4d34/r3 + r2d23/?'3) 
— (r4d34/r3 + r2d23/r3) 



Stage 
Global cost 
Local cost 



Cost 



y/nr2d2i/ri -^{d2i + d23) 



(^34 + r^dzA/ri) 



,13/14 



,19/14 



7WJT- 



Value 



Stage 



load 02001 load 02 o 03 load 03 o 04 



Global cost 
Local cost 



\J'r\ld2\ 



V^37^23 \/T3ld3A 



Cost 



Value 



7T077 



-WIT 



-T97T4- 



Stage 



load 05 



load 05 o 04 



Global cost 
Local cost 



\fnlrl 
\fnd^A 



Cost 
Value 



V?'4?'5 



5/2 



Vr2l-4rf2ici23rf,34 ^ ^ \/''2 ^21 ^231^34 ^ ^ 



,15/14 



□ 



The algorithms for finding fc-vertex subgraphs given in |13| [9] have complexity 0(n^ ''*) for finding a 
4-path, but it was not realized there that these algorithms apply to a much broader class of functions like 
associativity. The key property that is used for this application is that in the basic learning graph model 
the complexity depends only on the index sets of 1-certificates and not on the underlying alphabet. This 
property was previously observed by Mario Szegedy in the context of limitations of the basic learning graph 
model [8]. He observed that the basic learning graph complexity of the threshold-2 function is 0(n^/'^), 
rather than the true value 0(y^), as threshold-2 and element distinctness have the same 1-certificate index 
sets. 
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